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Abstract. In this paper, a regularity result for the Willmore flow is presented. 
It is established by means of a truncated translation technique in conjunction 
with the Implicit Function Theorem. 



1. Introduction 

The Willmore flow consists in looking for an oriented, closed, compact moving 
hypersurface T(t) immersed in R'^ evolving subject to the law 



\V{t) = -Ar(*)i?r(t) -2i?r(t)(i?^(,) -i^r(t)), 

lr(o) = ro. 



(1.1) 



Here V{t) denotes the velocity in the normal direction of F at time t and ^r{t) 
stands for the Laplace-Beltrami operator, while ^^r(t) is the normalized mean cur- 
vature of T{t). Finally, ^r(t) denotes the Gaussian curvature. 
The equilibria of (|l.ip appear as the critical points of the Willmore functional, or 
sometimes called the Willmore energy. For a smooth immersion / : F — > of a 
closed oriented two-dimensional manifold F, the Willmore functional is defined as: 

W{f)^ I Hj^^^da, (1.2) 
/(r) 

where da is the area element on /(F) with respect to the Euclidean metric in M.^. 
The critical surfaces of this functional, called the Willmore surfaces, satisfy the 
equation: 

^f{r)Hf(r) + 2-ff/(r) - 2-ff/(r)-f'^/(r) = 0. (1.3) 
The reader may consult [571 Section 7.4] for a brief historical account and a proof 
of this variational formula. The proof therein is derived by computing the critical 
points of all normal variations of the hypersurface /(F). 

A generalization of the Willmore functional (|1.2p in higher dimensions is studied 
by B.-Y. Chen [J. He extends (11.21) for smooth immersions / : F R^+i of the 
m-dimensional closed oriented manifold F into R''"+^: 



/(r) 
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with da standing for the volume element with respect to the Euclidean metric in 
Ijm+i^ The critical points of this functional are now of the form: 



m = 2, so this Euler-Lagrange equation agrees with (ll.Sp in the two-dimensional 
case. However, this generalization has the drawback that the corresponding Will- 
more functional is no longer conformally invariant except when m = 2. 

The Willmore problem has been studied by many authors, among them T.J. Will- 
more, W. Blaschke, B.-Y. Chen, J.L. Weiner, P. Li, S.-T. Yau, R. Bryant, R. Kus- 
ner, L. Simon, U.F. Mayer, G. Simonett, M. Bauer, E. Kuwert, R. Schatzle, 
U. PinkaU, I. Sterhng, M.U. Schmidt, CM. Fernando, and N. Andre. See for 



example [l[ll[i[I3[Il[Tl[Il[I3[T71[Il[Tl[2i[ll[Ml[13[2i[2Il. It is well- 



known that the Willmore functional is bounded below by An with equality only 
for the round sphere. Then the famous Willmore conjecture due to T.J. Willmore 
asserts that for any immersed 2-dimensional torus into M.^ we have W{f) > 27r^, 
and it suggests that the 2-dimensional Clifford torus achieves the minimum of the 
Willmore functional amongst all immersed tori in M"^. In 1982, P. Li and S.-T. Yau 
[T5] show that any immersion with W(f) < 8tt must in fact be an embedding. In 
other words, it will suffice to estimate W{f) for embeddings. A classification of all 
Willmore immersions / : S'^ — )■ IR.'^ is obtained by R.L. Bryant ^3j. The possible 



are Amr with n = 1, or n > 4 and n even, or ti > 9 and n odd. Existence and 
regularity for embedded tori in the Willmore conjecture has been proven by L. Si- 
mon [21], and later this result is generalized by M. Bauer, E. Kuwert |2] for an 
extension of the conjecture by R. Kusner [TT] to higher genus cases. An existence, 
uniqueness and regularity result on the Willmore flow is presented by G. Simonett 
[25] . It is proven therein that the Willmore flow admits a unique smooth solution. 
Moreover, this solution exists globally when it is initially close enough to spheres in 
the C^+"-topology and is exponentially attracted by spheres. In [T7], U.F Mayer 
and G. Simonett prove that the Willmore flow can drive embedded surfaces to a 
self-intersection in a finite time interval. Moreover, numerical simulations in |18| 
indicate that the Willmore flow can develop true singularities (topological changes) 
in finite time. E. Kuwert and R. Schatzle fT2l show that the smooth solutions are 
global as long as the initial Willmore energy is sufficiently small. Later, the same 
authors improve this result in |14j by finding an explicit optimal bound for the re- 
striction on the initial energy, that is, if the smooth immersion /g : F — >■ R'^ satisfies 
W{fo) < Stt, then the solution with initial data /o exists smoothly for all time and 
converges to a round sphere. Recently, in a breakthrough paper, CM. Fernando 
and N. Andre !¥| prove the Willmore conjecture for surfaces of arbitrary genus 
g > I, i.e., W{f) > Itv^ for all embedded F with genus <? > 1, and the equality 
holds iff F is conformal to the Clifford torus. 




values of 
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Assumptions: Throughout this paper, we always assume that {M,g) is a compact, 
closed, immersed, oriented, real analytic hypersurface in M.^ endowed with the Eu- 
clidean metric g with the exception of Section 3, wherein we remove the restriction 
on the dimension of M. The notation (-j-) always stands for the standard inner 
product in . We may find for M a normalized atlas (0^, (pK,)K.eA, where an atlas 
is said to be normalized if </Jk(Ok) = for all k E A. Here is the open unit ball 
centered at the origin in E^. Put tp,^ — f',^^- 

A family (7r„)„gA is called a localization system subordinate to (0^, (/5„)„gA if: 

(LI) ttk e Z^(Ok, [0, 1]) and (7r^)„gA is a partition of unity subordinate to the 
open cover (Ok)kgA- 

(L2) Any 7r„ and tt^ satisfying supp(7rK;) n supp(7r,,) ^ have their supports 
located within the same local chart. 

For any manifold satisfying the above assumptions, there exists a localization sys- 
tem. See \1, Lemma 3.2] for a proof. 

Notations: Throughout this paper. No stands for the set of all natural numbers 
including 0. For any time interval /, / always denotes the interior of /. 
Fix < a < 1. Let 7 e (0, 1] and Eq := /i"(M), Ei := h^+"{M). For notational 
brevity, we simply write 5^(0, M) and 5'(M,R) as ^(0) and 5^(M), where O is an 
open subset of and ^ stands for any of the function spaces in this paper. 

In the sequel, we always denote {Eo,Ei)y by E-y, where {■,-)-y is the continuous 
interpolation method. See [161 Definition 1.2.2] for a definition. Note that the 
continuous interpolation method (•, ■)g ^ coincides with (•, ■)g in the suggested ref- 
erence. In particular, we set {Eo,Ei)i := Ei. 

For some fixed interval / = [0, T] and some Banach space E, we define 

BUCi-M,E) := {uEC{i,E)-At^t^-'<u]e BUC{i,E), lim t^-'<\\u\\=Q}, 

\\u\\c,., := supii-'^h(t)||£, 
tei 

and 

BUC\_^{I,E) {u e C^ii,E) : u,ue BUCi-^{I,E)}. 

In particular, we put 

BUCail, E) := BUC{I, E) and BUCl{I, E) BUC^{I, E). 

In addition, if / = [0, T) is a half open interval, then 

Ci-,{I,E) := {v e C{i,E) : v G SC/Ci_-y([0, t], £:), t < T}, 

Cl_^{I,E) := {v e C\i,E) -.v^ve 

We equip these two spaces with the natural Frechet topology induced by the topol- 
ogy of BUCi-^{[0, t],E) and BUCl_^{[0, t],E), respectively. 
Last but not least, we set 

Eo(/) := C(/,£;o) and Ei{I) := C{I, Ei) D C^{I, Eq). 



4 



YUANZHEN SHAO 



It will be shown in this paper that the Willmore flow admits a real analytic 
solution jointly in time and space. Our motivation for a real analytic solution is 
mainly stimulated by the following facts: a compact closed real analytic manifold 
cannot have a " flat part" , and real analyticity in time implies that the hyersurface 
should move permanently in the interval of existence. 

Theorem 1.1. Let < a < 1. Suppose that Tq is a compact closed immersed 
oriented hypersurface in M.^ belonging to the class Then the Willmore flow 

has a unique local solution T = {T{t) : t £ [0,r)} for some T > 0. Moreover, 

M:= \J mxTit)) 

te(o,T) 

is a real analytic submanifold in R**. In particular, each manifold T{t) is real 
analytic for t £ (0, T). 

For any open subset O C M^, the little Holder space h'^(O) of order s > with 
s ^ N is the closure of BUC°°{0) in BUC'{0). Here BUC'iO) is the Banach space 
of all bounded and uniformly Holder continuous functions. The little Holder space 
h'^{M) on M is defined in terms of a smooth atlas, that is, a function u belongs to 
h%M) iff ^In^u e h''{M.^), for each k G A. 



2. Parameterization over a Reference Manifold 

In equation (|l.ip . if we fix an initial hypersurface Tq belonging to the class 
then by the discussion in [21] Section 4] we can find a real analytic compact closed 
embedded oriented hypersurface M, a function po e /i^+"(M) and a parameteriza- 
tion 

: M ^ R3, *p„(p) ■=p + pa{p)m{p) 

such that Tq = im(4'p„). Here I'Mip) denotes the unit normal with respect to a 
chosen orientation of M at p, and po • M — > {—a, a) is a real- valued function on M, 
where a is a sufficiently small positive number depending on the inner and outer 
ball condition of M. The reader may consult [HJ Section 4.1] for the precise bound 
of a. Thus To lies in the a-tubular neighborhood of M. In fact, it will suffice to 
assume Tq to be a C^-manifold for the existence of such a parameterization and a 
real analytic reference manifold. See [U] Section 4] for a detailed proof. 
Analogously, if T{t) is C^-close enough to M, then we can find a function p : 
[0, T) X M — > (—a, a) for some T > and a parameterization 

*p : [0,T) X M ^M3, ^p{t,p) ■.^p + p{t,p)m{p) 

such that r{t) — ini(5'p(i, •)) for every t G [0,r). It is worthwhile to mention 
that ^I'p admits an extension on R^, called Hanzawa transform, which was first 
introduced by E.I. Hanzawa in |10| . 

For any fixed t, I do not distinguish between p{t,-) and V'k(')) each local 
coordinate (0„,<^„) and abbreviate \['p(i, ■) to be Vfp := ^'p(t, •). In addition, the 
hypersurface r(t) will be simply written as Fp as long as the choice of t is of no 
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importance in the context, or p is independent of t. 
We put 

Here ||jo||^ := suppgy For any p G O, im(\I'p) constitutes a /i^+"-hypersurface 

Tp. In this case, defines a ft,^+"-diffeoniorphisni from M onto Tp. 

Here and in the foUowing, it is understood that the Einstein summation convention 
is employed and all the summations run from 1 to 2 for all repeated indices. 

In [21], J. Priiss and G. Simonett derive global expressions for many geometric 
objects of T p in terms of the function p. I will use some results therein to translate 
equation (|1.1[) into a differential equation in p. By [^H formula (23), (28)], we have 
the following explicit expressions for the components of the first fundamental form 
and the normal vector of T „: 



and 



9ij = 9t3 - "^Phj + P llhr + dipdjp, (2.1) 



yT= P{p){yM-a{p)). (2.2) 
In (|2.ip . the l^s are the components of the Weingarten tensor Lm of M with respect 
to g, i.e., Ly\ — Ifr"^ ^ tj, where {r^ = di} forms a basis of TpM at p S M and 
{t*} is the dual basis to {n}, i.e., (r'lrj) = 6^. The extension of Lm into R^, 
by identifying it to be zero in the normal direction, is denoted by L^, namely, 
— li'T'^ ^ Tj + ■ i^M ^E) i^M- It is a simple matter to check that 

t[ = {I - pLf^)n + md^p (2.3) 

forms the standard basis of T^^i^p^F^. In addition, the ^^^'s are the components of 
the second fundamental form L'^ of the metric g. Finally, gfj = (rf |rj) are the 
components of the first fundamental form of the Euclidean metric gr on Tp. We 
set G^{p) — {gfj)ij and G'p^(p) for its inverse. 

In (|2.2p . the terms a{p) and f3{p) read as 

a(p) = (/-pL^)-iVMP and /3(p) = [1 + |a(p)|2]-i/2. 
Here Vm is the surface gradient on M. 



For sufficiently small a > 0, (/ — pLfj^) is invertible. One can check that 



/ - pL^ = i^x - P^i^' ® + i^M ® I'M- 



Thus 



(I-pLf^) ^ =r{{py (g>Tj + iyM®i^M.. (2.4) 

where Rp = {rf{p))ij — [{Sj — plfjij]^^- By Cramer's rule, all the entries of Rp 
possess the expression 

'^'^ Qiip) 

in every local chart, where P/ and Ql are polynomials in p with real analytic 
coefficients and Ql ^ 0. 
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Substituting (/ - pL^^y^ by (|2^ . we get 

Hp)? = (r^(p)a,pr>i(p)5ipr'=) = g^^r^^{p)ri{p)d,pdip. 

Then 

m = [1 + |a(p)|2]-i/2 ^ [1 
Note that in every local chart 



where P^{p) is a polynomial in p with real analytic coefficients and Q^{p, djp) ^ 
is a polynomial in p and its first order derivatives with real analytic coefficients. 

The normal velocity can be expressed as 

V{t) = {dt^p\vv) = {ptvuW) = ?{p)pt- 
Therefore, the first line of equation is equivalent to 



Next we shall calculate the Gaussian curvature in terms of p. For simplicity, 

we write Kp instead of ^*pKYp ■ 

Because 

djn = ^ijTk + hjvu and djT'^ = -r^fc-r'' + l]vw, 
we may readily compute 
djLi, = a^/f Tfc - q^^f r' ® Tfe + r^.fc^f + ZjZf j^m ® + ® i^m • (2.5) 

Denote by = the second fundamental form of T p with respect to gr- Then 

by (|2.2p and (|2.3p . we can compute its components 1^^ as follows: 

= -((/ - pLl,)n + ^^M9.p|/3(9ji/M - d,a{p))) - {rl\^i^v) 

= p{kj + p{L^wT^\^Jm) + (T^|a,(VMp)) + ((/ - pLl^)T,\dj[{I - pL^)-1]Vmp) 
+ d,p{Md,[{I - pL^)-1]Vmp) + d,p{M{I - pLl,)-^[d,{Vup)])} 

= Pikj + pkk{T''\djVM) + (r,|9j(VMp)) + {n\dj{pLi^){I - pLl^y^Vup) 
+ dip{vw\dj{pLlj^){I - pLl^y^Vup) + dip{vu\dj{yup))} 

= P[hj - kki'^P + d,,p - T%dkP + ri(p)(9,zf + r^\/f - r^,/;;)p9,p 

+ ri{p)l':d,pdip + ri(p)Z^^/^p9,pazp + 
We have used (|2.5|) and the following facts in the above computation: 

• {m\dji'u) = 0. 

• (^fl^r) = 0. 
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. d,a{p) = {I- pLf^)-^d,{^up) + dj[{I - pL^)-i]VMp. 
. d,[{I - pL^)-i] = (/ - pL^)-ia,(pL^)(/ - pL^)-i. 
Therefore, dei(L^) can be expressed in every local chart as 

Here P^{p, djp, dijp) is a polynomial in p and its derivatives up to second order with 
real analytic coefBcients. Moreover, Q^{p) is a polynomial in p with real analytic 
coefficients. In particular, we have ^ 0. 

In full view of the above computations, within every local chart Kp ~ det[GY^{p)L^] 
can be expressed locally as 

_ ^ P^{p,djp,d,,p) 



As a straightforward conclusion of the above computation, we obtain an explicit 
expression for Hp := ^*iJrp: 

2Hp = gi^lf^ 

= (3{p)g]^ [h, - k.l^p + d.,p - T^^dkP + ri{p)l^djpdip 

+ r'kipKd.l'i + T%l1 - r':/^)pdip + ri{p)lpf^pd,pdip + p.pdup]. (2.7) 

The reader may also find a different global expression for Hp in ,21j formula (32)]. 
We can decompose Hp into Hp — Pi{p)p + Fi{p): 

Flip) = ^anh - hul'^p) - ^Tr[Gri(p)(L^ - pL^Lm)], 
where Tr(-) denotes the trace operator, and 

+ gWM'd.p + riipml + r%l1 - T':/,)p + riip)l';ilpd,pm 

in every local chart. Note that Tt[G^^{p)L^] changes like Hm under transition 
maps and thus is invariant. Analogously, we can check that Fi is a well-defined 
global operator. Hence so is Pi{p)- 

In addition, it is a well-known fact that \I'*Arp — Ap^**, where Ar^ and Ap 
are the Laplace-Beltrami operators on (Fpjgr) and (M,cr(p)), respectively. Here 
cr(p) := vE'p(7r stands for the pull-back metric of gr on M by 5'p. Then in every 
local chart, the Laplace-Beltrami operator Ap can be expressed as 

Ap^a'\p){d,dk-^]^{p)di). (2.8) 

Here (J^^{p) are the components of the induced metric (J*{p) of (7{p) on the cotangent 
bundle. Note that cr^*^ (p) involves the derivatives of p merely up to order one. 7*^ (p) 
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are the corresponding Christoffel symbols of cr(p), which contain the derivatives of 
p up to second order. 

There exists a global operator R{p) E C{h^^°' {M) , Eq) such that R{-) is well defined 
on 13 and: 



We set 



1 

1 . „ . . . 2 



Pip) - ^^pP^ + Pi-P)^ P e 



F{P) -^ApFi(p) + - J^Hp{Hl - K^), peUn h^+'^iM). 

Note that third order derivatives of p do not appear in F{p). Hence it is actu- 
ally well-defined on O. Based on the above discussion, these two maps enjoy the 
following smoothness properties: 

P eC°^{l3,C{Ei,Ei:,)) and e C°°(0, £:o). 

Definition 2.1. Let ^ G Nq. A linear operator A : I'(M) C(M) is called a 
differential operator of order I with continuous coefficients on M if for any u G V{M) 
it holds that 

for every local chart (0^, ^k) o,nd some differential operator Ak = 'l2\a\<i "^S^" with 
€ C(B^) defined on B^, and at least one of the Ak 's is of order I. In particular, 
when I = 0, Au = au for some a G C(M). 

By the above definition, P{p) is a fourth order differential operator with continuous 
coefficients on M for each p G 15. In every local chart {Of^,(pK,), the principal part 
of the local expression of P{p) can be written as 

P:ip) ■.^a'^\p)g'r'^^,kl. 
Given ^ G T*M, we estimate the symbol of P^{p) as follows. 

p:{pm-^'^*ipmo9ri^,o>c\e 

for some c > 0, and denotes the induced metric of gr on the cotangent bundle 
of M. Hence, P{p) is a uniformly elliptic fourth order operator acting on functions 
over M for each p e 13. By (23 Theorem 4.3], P{p) G n{Ei,Eo), namely that 
—P{p) generates an analytic semigroup on E^ with D{~P{p)) ~ Ei, p €l3. 



Now the Willmore flow (jl.ip can be rewritten as: 

ipt + P{p)p = Fip), 
\p{0)^Po, 

where po El3. A different characterization of the problem can be found in [51 [711^. 
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Applying O Theorem 4.1], the existence and regularity result in [25j can be restated 
as: 

Theorem 2.2. "25 > Theorem 1.1] Suppose that po e /i^+"(M). Then equation (f2J|) 
has a unique solution p such that 

P^c\ {J{po),Eo) n Ci (J{po),Ei) n c{J{po), h^+"iM)) n c^-z^^C J(po), e^,) 

2 2 

for any /3o S [0, ^]. Moreover, each hypersurface T{t) is of class C°° for t G J{p)- 



3. Parameter-Dependent DifFeomorphisms 

The main purpose of the last two sections is to show that the classical solution 
obtained in Theorem 12.21 is in fact real analytic jointly in time and space. To this 
end, I will construct a family of parameter-dependent diffeomorphisms acting on 
functions over M first. Because the construction applies to manifolds of arbitrary 
dimensions, in this section we assume that M is a m-dimensional manifold with the 
properties imposed in Section 1. 

For a given point p G M, we choose a normalized atlas (0^, fK)KeA for M such that 
ipiip) = G ffi.™. Choose several open subsets Bi in B'", the open unit ball centered 
at the origin in R™, in such a manner that: 

• Bi := B'"(0,ieo), for j = 1,2,3 and some eo > 0. 

• S3 CC Bi CCB'". 

Next, I further pick two cut-off functions on B™: 

• X G V{B2, [0, 1]) such that x\b^ = 1- We write Xk = flx- 

• C e V{Bi, [0, 1]) such that xl^., = 1- We write Ck = '^*C- 

We define a rescaled translation on B™ for any p, G B(0, r) C M™ with r sufficiently 
small: 

This localization technique in Euclidean spaces was first introduced in [S' by J. Es- 
cher, J. Priiss and G. Simonett to establish regularity for solutions to parabolic and 
elliptic equations. 

Given a function v G Li^;oc(B'"), its pull-back and push-forward induced by 9^^ are 
defined as: 

eiv:=voe^, and 9^:^ v o Q-"' . 
The diffeomorphism 0^ induces a transformation on M by: 



q qiOi. 
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It can be shown that 6^ G Diff °°(M) for /i e B(0, r) with sufficiently small r > 0. 
See [23] for details. 

For any u e ii.ioc(M), we can define its pull-back and push-forward induced by 0^ 
analogously as: 

e*M:='uoep and O^u := u o O''^ . 
We may find an explicit global expression for the transformation 0* on M, 

e;u ^ ^l9;ri{Ciu) + {1 - Ci)^. 

Here and in the following it is understood that a partially defined and compactly 
supported function is automatically extended over the whole base manifold by iden- 
tifying it to be zero outside its original domain. 

Likewise, we can express 0» as 

ei: ^ ipi9!:riiCiu) + {1 - Ci)u. 

Let / = [0, T], T > 0. Assuming that J C (0, T) is an open interval and to G J is a 
fixed point, we choose eo to be so small that M{tQ, 3eo) C J. Next we pick another 
auxiliary function 

eeP(B(<o,2eo),[0,l]) with C|B(to,so) = 1- 

The above construction now engenders a parameter-dependent transformation in 
terms of the time variable: 

g^{t) := t + ^(i)A, for any t e / and A G M. 

Now we are in a situation to define a family of parameter-dependent transformations 
on / X M. Given a function u : / x M — > M, we set 

where T^(i) = e*^^)^ and (A,^) G B(0,r). 

It is important to note that ^^,^^(0, •) — m(0, •) for any (A,/i) G B(0,r) and any 
function u. 

The importance of this family of parameter-dependent diffeomorphisms lies in the 
following theorems. Their proofs as well as additional properties of this technique 
can be found in [53]. 

Theorem 3.1. Let fc G Nq U {oo, uj}. Suppose that u G C{I x M). Then we have 
that u G C^{I X M) iff for any {to,p) G / x M, there exists r — r{to,p) > and a 
corresponding family of parameter- dependent diffeomorphisms 0^ ^ such that 

[(A,^) ^ei^ii] GC'^(B(0,r),C(/x M)). 
Here to is the symbol for real analyticity. 
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Proposition 3.2. Suppose that u E Ei(/). Then ux.fj. G lEi(^), and 
dt[uxj = (1 + C'A)ei,^ut + Ba.mK,^), 

where 

Furthermore, B\ o — 0. 

Proposition 3.3. Let s £ [0, /i] and I e Nq. Suppose that A is a differential 
operator of order I with continuous coefficients on M satisfying a'^ G _BC''(B™) and 
al e SC''(B™) nC"(0) for some open subset such that CC CC B". Then 

[fi^T^AT-^] e C-(B(0,r),C(/,/:(/i^+'(M),/i^(M)))). 

Proposition 3.4. Let s > 0. Suppose that u £ C"(-(/'i(0)) n h''{M), where is 
defined in Proposition \3.3i Then 

[fi^T^u]eC'-{M{0,r),C{L,h%M))). 



4. Real Analyticity 



By setting G{p) := P{p)p — F{p), we may rewrite equation p.9p as 

1p(0) = Po. ^'^ 



Theorem 4.1. Let < a < 1. Suppose that pa e /i^+"(M). T/ien equation (|4l" 
/las a unique local solution p in the interval of maximal existence J{po) such that 

pe C"(j(po) X M). 



Proof. I will indicate herein all the key steps of the proof. More details can be 
found in [53]. 

For any {to,p) G J(po) x M and sufficiently small r > 0, a family of parameter- 
dependent diffeomorphisms 9^^ can be defined for (A,/z) S B(0,r). Henceforth, 
we always use the notation p exclusively for the solution to (|2.9I) and hence to (|4.1I) . 
Set u := px^ii. Then as a consequence of Proposition 13.21 u satisfies the equation 

= -il + emi^G{p)+BxAu) 
= -il + e\)T^Giglp) + BxAu) 
= -(1 + rA)T^G(r-iu) + BxA^) -Hxju). 
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Pick / : [£,r] CC J{pq) such that to E I and B(to,3£o) CC /. Then we define 
Eo(/) and as in Section 1 by moving the initial point from to e. Set 

E?(/) ■.= {veEi{I) : \\v\\oo<a}, 
where Halloo := sup(t 9)1- 

For ^ G H{Ei, Eq), we say that (Eo(/), Ei (/)) is a pair of maximal regularity of 
A, if 

(^+A7e) eIsom(Ei(/),Eo(/) xi;i), 
where 7^ is the evaluation map at e, i.e., jsiu) — u{e). Next we define 

$ : E«(/) X B(0,r) ^ Eo(/) x E^ as a>(z;, (A,/.)) ^ ht^^'']""} 

Note that = Q for any {X, fi) e B(0,r). 

(i) My first goal is to prove that $ e C"(E!f(/) x B(0,r),Eo(/) x Ei). 

By Proposition O Ba,m e C"(B(0, r), C(/, £(£'1, i;o))). We define a bilinear and 
continuous map: 

/ : C{I,CiEi,Eo)) X El (J) ^ Eo(/), (T{t),u{t)) ^ T{t)iuit)). 

Hence [{v,{X,fi)) ^ f{Bx,^.,v) = Bx,^{v)] e C"(EJ(/) x B(0, r), Eo(/)). 
On the other hand, let vr = J2rjec{i) '^v'' ''^here 

C(f ) {77 e A : supp(7r^) n supp(7ri) ^ 0}. 
We decompose G into 

G = ttG + I],,^c(i) 

According to our construction of 0* and of the localization system, we may assume 
that 7r|o = 1, where is defined in Proposition [331 with m — 2. See [1] Lemma 3.2] 
for details. 

Taking into account ()2.6p . p.7p and p.Sp . in every local chart (0^, (^k) and for any 

V e Ei(/), G(w) can be expressed as 

/3^^(t;)pg(t;,--- ,d,jkiv) 
det{G^{v)Y^dei{[a{v)]Y^QG{vy 

where h,si,S2 S N. [ct{v)] is the matrix representation of the metric (t{v). Here 
cr(v) is defined in a similar manner to cr(p) with p replaced by v. Analogously, 
{v) is defined in a similar way to G^{p). Meanwhile, is a polynomial in 

V and its derivatives up to fourth order with real analytic coefficients, and Q*^ is 
a polynomial in v with real analytic coefficients. In particular, det([cr(t;)]) only 
involves first order derivatives of v. 
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Therefore, TrG{v) can be decomposed globally into 

Vo + rlv-Vl^v + --- + V{v ■ ■ ■ riv 
Qo + Q\v---Qlv+--- + Qlv ■ ■ ■ Qfv ' 

where 7^0, Qo S C°°(M)nC"(?/'i(0)). The T'j/s are Imear differential operators with 
continuous coefficients on M up to fourth order, and the Qj. 's are linear differential 
operators of order at most one with continuous coefficients on M . Their coefficients 
in every local chart satisfy that ag e BC°°{B^) and £ BC°°{B^) n C"^(0). By 
Proposition we deduce that 

{T^r„,T^Qo)] eC^{M{0,r),C{I,Ei) X C{I,Ei)). 

Analogously, it follows from Proposition 13.31 that 

[/i ^ T^riT-'] G C"(B(0,r),C(/,/:(i?i,i?o))) 

and 

ifi^T^QlT-'] e C"(B(0,r),C(/,/:(i?i,/z3+"(M)))). 

Combining the above discussion with point-wise multiplication theorems on Rie- 
mannian manifolds, we infer that 

[{v,fi) ^ T^{7rG)T-'v] e C^iEUl) x B(0, r), Eo(/)). 

Applying these arguments repeatedly to the other terms tt^G, we conclude that 

$ e C"^(E?(/) X B(0,r),Eo(/) x Ei). 



(ii) Next we look at the Frechet derivative of $ in the first component 
(A, 



wt + il+ i'X)T^DG{T-^v)T-^w - Bx.^,{w) 

Thus 

Observe that DG{p) is a fourth order linear differential operator whose coefficients 
satisfy a'^ & Eq. The principal part of DG{p) in every local chart coincides with that 
of P{p), that is, P^{p). By the discussion in Section 2, we know that DG{p{t, •)) 
is a uniformly elliptic operator for every fixed t > 0. As a consequence of |23[ 
Theorem 4.5, Proposition 4.7], it follows that (Eo(/),Ei(/)) is a pair of maximal 
regularity for DG{p{t, •)). 
We set Ait) = DG{p{t, ■)). It follows that 

+A{s),je) e Isom(Ei(/),Eo(/) X Ei), for every s e I. 
By O Lemma 2.8(a)], we have 

(-^+^(-),7e) eIsom(Ei(/),Eo(/) xi;i). 
at 

Now we are in a position to apply the Implicit Function Theorem. It follows right 
away that there exists an open neighborhood, say B(0, ro) C B(0, r), such that 

[(A,m)h^Pa,m] eC"(IB(0,ro),Ei(/)). 
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As a consequence of Theorem 13.11 we deduce that p G C'^{J{po) x M). This com- 
pletes the proof. □ 

Proof of Theorem \l.l\ For each (to,q) E M = Ute,7(po)('f^J' ^ ^(^))' there exists a 
p e M such that *p(to,p) = (?• Here T{t) = im(*p(t, •)). Theorem BTTl states that 
there exists a local patch (0^, v^k) such that p G 0^ and p o ip^^ is real analytic in 
J{po) X B^. Therefore, we conclude that 

This proves the assertion of Theorem 11.11 □ 
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